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I. INTRODUCTION 

The spectrum of known cs states can be described as 
two S-wavc states (-Dj", -D*+) with = 0^, 1^, and four 
P-wave states (15*0(2317)+, 15,1(2460)+, 15,1(2536) + , 
L>;2(2573)+) with = 0+,l+,l+,2+. Whether this 
picture is correct remains controversial because the states 
at 2317 McV/c^ and 2460 MeV/c^ [1] had been expected 
to lie at much higher masses [2]. Since the discovery 
of the new Df mesons much theoretical work has been 
done, however new experimental results are needed in 
order to understand this sector of spectroscopy. 

Recently, two new _D+ mesons have been discovered, 
D*i(2710j+ [3, 4] and i:>*;(2860)+ [3]. The analysis of 
D*i(2710)+ produced in B decays gives the assignment 
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JP = 1-. For D*j(2860)+ assignments of = 0+ [5, 6] 
and jP = 3~ [7] have boon proposed. 

We report here on a search for new D^j mesons in 
the mass spectrum of D'-*^K inclusively produced at 
the PEP-II asymmetric-energy e~^e~ storage rings and 
recorded by the BaBAR detector. This paper is organized 
as follows. In Sect. II we give a short description of 
the BaBar experiment and in Sect. Ill we describe the 
data selection. Section IV is devoted to the study of the 
DK system and in Sect. V we present the study of the 
D*K system. In Sect. VI we describe fits to the D*K 
mass spectrum while in Sect. VII we present an analysis 
of the angular distributions. Measurements of ratios of 
branching fractions are described in Sect. VIII and we 
summarize the results in Sect. IX. 
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FIG. 1: (a) K tt"*", (b) K Tr+Tr"*" mass distributions for 
candidate events in channels (1) and (2) respectively. 
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II. THE BABAR EXPERIMENT 

This analysis is based on a 470 fb~^ data sample 

recorded at the T(4S') resonance and 40 MeV below the 
resonance. The BaBAR detector is described in detail else- 
where [8] . We mention here only the parts of the detector 
which are used in the present analysis. Charged particles 
are detected and their momenta measured with a combi- 
nation of a cylindrical drift chamber (DCH) and a silicon 
vertex tracker (SVT), both operating within the 1.5 T 
magnetic field of a superconducting solenoid. Informa- 
tion from a ring-imaging Cherenkov detector combined 
with energy-loss measurements in the SVT and DCH 
provide identification of charged kaon and pion candi- 
dates. The energies and locations of showers associated 
with photons are measured with a CsI(Tl) electromag- 
netic calorimeter. 



III. DATA SELECTION 

We reconstruct the inclusive processes [9] listed in Ta- 
ble I. A particle identification algorithm is applied to all 
the tracks. Charged kaon identification has an average ef- 
ficiency of 90% within the acceptance of the detector and 
an average pion-to-kaon misidentification probability of 
1.5% per particle. 

For all channels we perform a vertex fit for the D*^ and 
daughters and require a probability greater than 
0.1%. For the tt'^ candidates in channels (3), (4), and 
(6), we combine all photons with energy greater than 
30 MeV in pairs, perform a fit with a tt" mass con- 
straint, and require a probability greater than 1%. 
For the K"n^Ti^ decay channel we also perform 

a kinematic fit with a £>° mass constraint. We obtain 
K'g — > 7r+7r~ candidates by means of a vertex fit and re- 
quire a probability greater than 2%. We accept only 
candidates with decay length greater than 0.5 cm. To 
obtain Z)*+/f° candidates, where — > £)°7r+, we com- 
bine fitted and candidates with a tt"*" candidate 
using a vertex fit which constrains the overall vertex to be 



located in the interaction region, requiring a probabil- 
ity greater than 0.1%. Similarly, for D*^K^ candidates 
where D*+ D+tt° we combine fitted D+, K° and tt° 
candidates using a vertex fit which constrains the overall 
vertex to be located in the interaction region, requiring 
a probability greater than 0.1%. Background from 
e+e" BB events is removed by requiring the center 
of mass momentum p* of the DK or D*K system to be 
greater than 3.3 GeV/c. 

To improve the signal to background ratio for channels 
with D° K~TT^ we study the distribution of the an- 
gle 6k- formed by the K~ from D° decay in the K~Tr~^ 
rest frame with respect to the K~Tr~^ direction in the 
laboratory system. This distribution is expected to be 
fiat. We observe an accumulation of combinatorial back- 
ground close to cos6x- = 1- We improve the signal to 
background ratio by requiring cos^^- < 0-9- 

To improve the signal to background ratio for £)+ —> 
K~7:^t:^ , we compare the three-momentum and its 
fiight direction and define d^y as the signed projected 
distance in the transverse plane. Background events are 
removed by requiring d^y > 0. The resulting i^~7r+ and 
i^T" Tr+Tr"*" invariant mass spectra for candidates in chan- 
nels (1) (where we require in addition a reconstructed 
and (2) (where we require in addition a recon- 
structed Kg) are shown in Fig. 1. There are on average 
1.01 candidates per selected event in both samples, and 
all candidates are retained for further analysis. 

We fit the K~n^ and K~'K~^n~^ invariant mass spec- 
tra using a linear background and a single Gaussian peak 
obtaining ctuo = 7.6 MeV/c^ and aD+ = 6 MeV/c^. The 
signal region is defined within ±2cr while sideband re- 
gions are defined within (— 6cr, — 4cr) and (4(j, 6(t). The 
signal region contains 1.98 xlO^ combinations with 
a purity P = Ns/{Ns + Nb) = 0.84, where Ns (Nb) 
is the number of signal (background) combinations. The 
D"*" signal region contains 0.58 xlO^ combinations with 
a purity P = 0.75. 
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TABLE I: List of reconstructed final states. Here X indicates the rest of the event, with any number of charged or neutral 
particles. 



Channel Reaction D* decay mode D decay mode 



(1) 


e^ 


"e^ 


D^K+X 






^ R- 






(2) 


e^ 


~e~ 


D+KlX 




D+ 


K- 




%+ 


(3) 


e^ 


~e~ 


_^ D*°K+X 




D» 


R- 






(4) 


e^ 


~e~ 


D*+KlX 


D*+ ^ Z)+7r° 


D+ 


R- 


~n' 


%+ 


(5) 


e^ 


~e~ 


W+KlX 


D*+ ^ 13%+ 


D° 


R- 






(6) 


e^ 


~e~ 


D'+KlX 


D*+ ^ D"7r+ 


73" 


R- 






(7) 


e^ 


~e~ 


D*+KlX 


D*+ ^ £)%+ 




R- 




- + — 



IV. STUDY OF THE DK SYSTEMS 

We first study the D^K^ and D^K^ mass spectra. In 
an inclusive environment the and £>+ can come from 
D* decays. Candidate D^K+ pairs where the is a 
Z3*-dccay product arc identified by forming Z3"7r+, L'^tt'^, 
and combinations and requiring that the invariant- 
mass diff'erence between one of those combinations and 
the be within ±2(7 of the known D* — D mass differ- 
ence. Events belonging to these possible reflections (ex- 
cept for D*^ —> D^j events, which could not be isolated 
cleanly) are removed. In the same way, D'^K'^ combina- 
tions where the £>+7r'' and Z)+ invariant-mass difference 
is found to be within ±20- of the known D* — D mass 
diflference axe removed. 

We also study the distribution of 6ji+ {(^k% ) i the angle 
between the K'^ {K^^) direction in the DK rest frame 
and the DK direction in the laboratory frame. We ex- 
pect the distribution of this angle to be symmetric around 
zero [10], but we observe an accumulation of combinato- 
rial background close to cos 9k = — 1- Due to the jetlike 
nature of the reaction e+e" — > cc we interpret this back- 
ground as due to combinations for which the K comes 
from the jet opposite to the D meson. We therefore ap- 
ply a conservative cut requiring cosOk > —0.8. 

The resulting D^K^ and D+if" mass spectra are 
shown in Fig. 2. To improve the mass resolution, the 
nominal D mass and the reconstructed 3-momentum are 
used to calculate the D energy for channels (1) and 
(2). The two mass spectra in Fig. 2 present similar fea- 
tures. The single bin peak at 2.4 GeV/c^ results from 
decays of L'si(2536)+ to D*^K+ or D*+Kl in which 
the 77° or 7 from the D* decay is missed. Since the 
£'si(2536)+ is believed to have = 1+, decay to DK 
is forbidden by angular momentum and parity conserva- 
tion. We also observe a prominent narrow signal due to 
the 13*2(2573)+, a broad structure centered at the mass 
of the Z3*i(2710)+, and a narrower structure at the posi- 
tion of the £)*j(2860)+. 

We perform a simultaneous binned to the 

two sideband-subtracted DK mass spectra shown in 
Fig. 3(a) and Fig. 3(c). The fit range extends from 
2.42 GeV/c^ to 3.2 GcV/c^ (the lower bomid is chosen 
to exclude the Dsi(2536)+ reflection). The background 



for each of the two DK mass distributions is described 
by a threshold function: [m — mth)" e"''™"^™ 
where rrith = mo -\- rriK- In this fit, the 13*2(2573)+, 
£»;i(2710)+ and Z3;j(2860)+ peaks are described with 
relativistic Breit-Wigner lineshapes where spin-2 is as- 
sumed for 23*2(2573)+, spin-1 for £)*i(2710)+ and spin- 
for D*j{2860)+. The Breit-Wigner function includes 
Blatt-Weisskopf form factors [11]. Mass resolution es- 
timates are obtained using simulated events. We obtain 
single-Gaussian ct's of 2.7 MeV/c^ and 3.6 MeV/c^ at DK 
masses of 2.71 and 2.86 GeV/c^ respectively. Since the 
width values for the resonances present in the DK mass 
spectra are much larger than these, resolution eff'ects are 
ignored. 

The result of the fit is shown in Figs. 3(a) and 3(c), 
and the parameter values obtained arc summarized in Ta- 
ble II (Fit A). Figures 3(b) and 3(d) show also the fitted- 
background-subtracted D'^_^+K~^ , and in- 
variant mass distributions. The fitted parameter values 
for the D* J (2860)+ state are in agreement with our pre- 
vious measurement [3], while the central value of the 
23*1(2710)+ mass is slightly higher than before. 



V. STUDY OF THE D*K SYSTEM 

The Am — m(Z37r) — m,{D) distributions for the five 
channels (3)-(7), are shown in Fig. 4. Backgrounds are 
small for channels (5)-(7) but larger for channels (3)-(4). 
Table III gives the fitted parameter values of the Am 
distributions together with purities and the definitions 
of signal and sideband regions. The values of a reported 
in Table III are obtained from fits performed using a 
polynomial background and a single Gaussian. 

We have also studied the distributions of the angle 
{6*j^o ) between the K+ {Kf) direction in the D*K rest 
frame and the D*K direction in the laboratory frame. 
We expect the distributions of this angle to be symmetric 
around zero for signal but we observe an accumulation 
of combinatorial background close to cos 6*j^ = — 1. As 
in the case of the DK system, we interpret this as being 
due to combinations for which the K comes from the jet 
opposite to that yielding the D* meson. We therefore 
apply the conservative selection criterion cos0|f > —0.8. 
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FIG. 2: DK invariant mass distributions for (a) D'^_^^K^ , and (b) D'^_ ^^^_f_K^^ . Shaded histograms represent the D mass 
sideband regions. The arrows indicate the expected positions of the 1)^2(2573)"*', Z)Ji(2710)"'", and -D*j (2860)"^ . 
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FIG. 3: Sideband-subtracted Dif invariant mass distributions for (a) D^_^^K^ , (c) _Dj_^^^^/('s ; (b) and (d) show the 
fitted-background-subtracted mass spectra. The curves show the functions described in the text. 
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TABLE II: The xVNDF and resonance parameter values obtained from the fits to the DK and D* K mass spectra. Masses 
and widths are given in units of MeV/c^ and MeV, respectively. Uncertainties are statistical only. 



Fit 


xVndf 




D,* J (2860)+ 




(3040) + 


A [DK) 


85/56 


m = 2710.0 ± 3.3 


m = 2860.0 ± 2.3 










r = 178 ± 19 


F = 53 ± 6 






B {D*K) 


51/33 


m = 2712 ± 3 


m = 2865.2 ± 3.5 


m = 


3042 ± 9 


r = 103 ± 8 


F = 44 ± 8.3 


F = 


214 ± 34 


C {DK + D*K) 


147/91 


m = 2710 ± 3 


m = 2860 ± 2 


m = 


3045 ± 8 






F = 152 ± 7 


F = 52 ± 5 


F = 


246 ± 31 








m = 2866 ± 3 












F = 43 ± 6 






D [DK + D'K) 


149/93 


m = 2710 ± 2 


m = 2862 ± 2 


m = 


3044 ± 8 






F = 152 ± 7 


F = 48 ± 3 


F = 


239 ± 35 


E (D*K) 


65/38 


m = 2716.7 ±2.5 












F = 108 ± 5 








F {D*K) 


39/34 






rn = 


3047 ± 12 










F = 


2i() ± 46 




0.12 , ,0.18 

Am GeV/c' 
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FIG. 4: The Am distributions for channels (3)-(7) after applying the corresponding D-candidate mass selection criteria. 
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TABLE III: Fitted parameters of the Am distributions together with purities and definitions of the regions used for signal and 
background. 



Channel 


mass a 
MeV/c^ MeV/c^ 


purity (%) Signal region 


Sideband 


region 


(3) Am(Z?^_^+^«) 


142.02 


1.08 


83.3 


±2.5(7 


10(7 - 


15a 


(4) Am(0+_^+^+^«) 


140.63 


0.893 


76.6 


±2.5ct 


10(7 - 


15ct 


(5) Am(D^_^+^+) 


145.43 


0.288 


94.9 


±5(T 


12a - 


22ct 


(6) Am{D%^^^^„7v+) 


145.43 


0.351 


87.1 


±3cr 


10(7 - 


16ct 


(7) Am(7?^_^+^+^_^^ 


-) 145.43 


0.266 


90.5 


±5(7 


12(7 - 


22a 
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2.5 3 

m(D*K) GeV/c' 
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FIG. 5: The m{D*K) distributions for channels (3)-(7); the shaded histograms show the mass spectra from the D* Am sideband 
regions. The arrows indicate the expected positions of the Dsi(2710)"'" and D*j(2860)'^. 



Figure 5 shows the reconstructed D*K mass spectra 
for channels (3)-(7) of Table I. The shaded distributions 
represent the background estimates from the Am side- 
band regions. Each distribution shows a narrow spike at 
threshold due to the Z)si(2536)"'" meson. We also observe 
structures around 2.71 and 2.86 GeV/c^. 

We have compared these mass spectra with those from 
generic e~^e~ — > cc Monte Carlo events. These events 



were generated using a detailed detector simulation and 
subjected to the same reconstruction and event-selection 
procedure as used for the data. We find that the sim- 
ulation underestimates the Dsi(2536)"'" and 15*2(2573)"'' 
signals relative to the background. No such discrepancy 
is found in the study of non-strange final states, therefore 
we attribute this effect to poor knowledge of the strange- 
charmed meson cross sections [12]. We apply weights to 
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the Dsi{2536)+ and D*2(2573)+ production in the Monte 
Carlo in order to obtain better agreement with the data. 

Figure 5(3) shows the presence of a peaking back- 
ground in channel (3) around 2.7 GeV/c^. Using the 
Monte Carlo data we identify this reflection, which is 
present in the signal and the sideband regions, as being 
due to the £>*2(2573)+. Combinations of DK originat- 
ing from this narrow peak associate with a random 7r° to 
produce a relatively narrow structure in the 2.7 GeV/c^ 
region. Our Monte Carlo study verifies that this reflec- 
tion is almost completely removed by performing the Am 
sideband subtraction. 

The total D*K mass spectrum, Aw-sideband- 
subtracted and summed over channels (3)-(7), is shown 
in Fig. 6 and compared with that obtained from Monte 
Carlo simulations. The D*K mass spectrum, above 
the £'si(2536)+ signal, shows the presence of structures 
around 2.71, 2.86, and 3.04 GeV/c^. Corresponding reso- 
nance contributions are not included in the Monte Carlo 
simulation. Since such enhancements are not evident in 
the Monte Carlo D*K mass spectrum, we conclude that 
these structures are not produced by reflections from 
known resonances. Monte Carlo simulations also show 
that these enhancements are not due to reflections from 
the Dsj resonances observed in the DK mass spectrum. 

A structure close to 2.57 GeV/c^ is seen in the Monte 
Carlo and is due to the decay Dl^{2b72>)+ D*K, in- 
cluded in the simulation. However the data do not show 
evidence for such a decay. 



those of the £>*j(2860)+ obtained from the DK mass 

spectrum. 

We next repeat the fits removing the resonances one 
by one from the PDF. We compute the statistical sig- 
nificance of each structure as yAx^ where A^^ is the 
difference in the fit with and without the resonance 
included, taking into account the variation in the number 
of parameters (ANDF = 3). We obtain statistical signif- 
icances of 12.4, 6.4, and 6.0 standard deviations for the 
£>;i(2710)+, £>*j(2860)+, and £)^j(3040)+ respectively. 

We then perform simultaneous fits to the two DK mass 
spectra and to the total D*K mass spectrum to better 
constrain the D*i(2710)+ parameters. We first test the 
possibility that the structure around 2.86 GeV/c^ in the 
D*K mass spectrum is different from the D*j(2860)+ 
state observed in the DK mass spectrum by adding to 
the fit two new parameters. The results from the fit are 
summarized in Table II (Fit C). We find the parameters 
of the Z?*j(2860)"'' in the DK mass spectrum consistent 
with those measured in the D* K mass spectrum. 

Assuming therefore that we are observing the same 
state, we constrain the £'*j(2860)"'" parameters to be the 
same in both the DK and D*K mass spectra (Fit D). 
Taking this as our reference fit, we obtain the following 
parameters for the three states: 

m(i?:i(2710)+) = 2710 ± 2stat(;7')syst MeV/c^ 
r=149±7stat(ti2)syst MeV, 



VI. FITS TO THE D' K MASS SPECTRUM 

We perform a binned minimum fit to the combined 
D*K mass spectrum. The fit is performed in the region 
(2.58-3.48) GeV/c^. The background is parametrized as 
g-/3im-7m -5m ^ ^\^{(^ provides a good description of the 
Monte Carlo in the same mass range. The D^j peaks are 
described with relativistic Breit-Wigner lineshapes. Here 
and in the following fits we assume ~ \^ and — 
respectively for D;i(2710)+ and D;j(2860)+ decays to 
DK. For the D*K system, we use an angular momentum 
L = 1, L = 3, and L = for D*i(2710) + , L'*;(2860)+, 
and £'sj(3040)"'" respectively. Average mass resolutions 
are 2.5 MeV/c^ and 3.5 MeV/c^ at D*K masses of 2.71 
and 2.86 GeV/c^ respectively. Since the width values 
for the resonances present in the D* K mass spectra are 
much larger than these, resolution effects are ignored. 
We observe the presence, above the £'*i(2710)"'" and 
J (2860)"*", of a new broad structure peaking at 3.04 
GeV/c^. The resonance parameters resulting from the fit 
are given in Table II (Fit B) and the corresponding fitted 
curves are shown in Fig. 7. Modifying the background to 
include an extra term in the exponential docs not im- 
prove the fit significantly. We find that the width of the 
£)*i(2710)+ differs somewhat between the DK and D*K 
fits, while the parameter values for the structure at 2.86 
GeV/c^ in the D*K mass spectrum are consistent with 



m(D:^(2860)+) = 2862±2,tat(l2)syst MeV/c^ 

± ^stat ± 6,yst MeV, 
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m(D«j(3040)) = 3044 ± 8stat(lf )syst MeV/c^ 
r = 239±35stat(ll2)syst MeV. 



(2) 



(3) 



Here systematic uncertainties take into account the 
range of values obtained in different fits, including fits 
to the spectra obtained after modifying the p* selection 
criterion, modifying the Arn. criteria, and removing the 
cos 9k requirement. They also account for uncertainties 
in the spin assignment. 



VII. ANGULAR ANALYSIS 

Since we observe both Dl^{2710)+ and i:>*^(2860)+ de- 
cays to both DK and D*K, we assume they have natural 
parity = 1~, 2+, 3~, ... (J^ = 0+ is ruled out because 
of the D*K decay). We further test this hypothesis us- 
ing angular analysis. We compute the helicity angle 9^ 
as the angle formed by the n from the D* decay with 
respect to the kaon, in the D* rest frame. The angular 
distribution for natural parity is expected to be [13] 



dN 



d cos 9h 



^ 1 - cos^ 9h , 



(4) 



12 



4000 

O 

> 3000 



S 2000 



I 1000 







V 



(a' 



2.5 3 

m(D'K) GeV/c' 




2.5 3 

m(D'K) GeV/c' 



FIG. 6: The sideband-subtracted m[D* K) distributions sunmiod over channels (3)-(7). (a) Data, (b) Monte Carlo. The arrow 
near threshold indicates the position of the peak due to the Z)si(2536)"'", which is off-scale. The other arrows indicate the 
expected positions of the Z)*i(2710)+, D*j(2860)+, and Dsj(3040)+. 
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FIG. 7: (a) Fit to the D*K invariant mass spectrum, (b) Residuals after subtraction of the fitted background. The curves are 
described in the text. 



since, for the parity and angular momentum conserving 
decay of such a parent state, the coupling in the par- 
ent rest frame to the helicity-0 D* state involves a van- 
ishing Clebsch-Gordan coefficient. Figure 8 shows the 
D*K mass spectrum separated for |cos6';j| < 0.4 ((a),(b)) 
and |cos6';i| > 0.4 ((c), (d)). We clearly observe an en- 
hanced signal to background rate for the _D*]^(2710)+ in 
the |cos^?;i| < 0.4 region. This does not hold for the 
£>^j(3040)+. The non-observation of L»sj(3040)+ DK 
also suggests unnatural parity for this state. 

The mass spectra separated according to the value of 
\cos9h\ allow a better determination of the Dl-^{271Q)^ 



and Ds,7(3040)+ parameters. When fitting the |cos0/i| < 
0.4 data (Fit E) wc fix the D*j(2860)+ and L>,j(3040)+ 
shape parameters to those resulting from the simultane- 
ous fit of the DK and D*K mass spectra (Fit D). In 
fitting the jcos^^j > 0.4 data (Fit F), we fix the param- 
eters of 15*1(2710)+ and 15*^(2860)+ to those from Fit 
D. The resulting D*i(2710j+ and i:>,,,/(3040)+ parame- 
ters are given in Table II and the fit results are shown by 
the curves in Fig. 8. 

We have studied the cos^?t distributions for the 

r>|i(2710)+ and D*j(2860)+ by producing D*K mass 
spectra in six intervals of cos Oh- The mass spectrum in 
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FIG. 8: Fits to the D*K invariant mass spectra for (a) |cos0/i| < 0.4; (c) |cos^h| > 0.4. (b) and (d) show the data after the 
fitted background is subtracted. 



each interval was fitted with fixed resonance parameters. 
However these values have all been varied within their 
statistical and systematic errors. The efficiencies as a 
function of cos 9h in the two mass regions are obtained 
from Monte Carlo simulation of the five channels involved 
in the analysis. Wc find that the efficiency is almost uni- 
form as a function of cos 9h with a small slope which we 
parametrize by a linear function. 

The efficiency corrected i:»;i(2710)+ and D;j{2860)+ 
yields are plotted in Fig. 9, together with the normalized 
expectations for natural parity. The curves have x^/NDF 
of 18.7/5 and 6.3/5 respectively. The large obtained 
for the D*^{2710)~^ is related to the large uncertainties in 
the background parametrization. Other spin hypotheses 
have been tested but they give much larger values. 
We conclude that both states are consistent with hav- 
ing natural parity. We do not perform a similar analysis 
for £>sj(3040)+ because of the large uncertainties arising 
from fitting a very broad resonance with limited statis- 
tics. 



3 I 1 I 




-1 1-1 

cos -tfi, COS 

FIG. 9: Distributions in cosOh for (a) the i)*i(2710)+ and (b) 
the 7?: J (2860)+. 
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VIII. BRANCHING FRACTIONS 



Prom Table I it can be seen that it is possible to ob- 
tain ratios of branching fractions with reduced systematic 
uncertainties, for D*^{2710)+ and D*j{2860)+ by using 
channels (3),(1) and (4), (2) respectively. 

These ratios are computed as: 



r-4 



NjPlj D*K) e{D*j DK) 
N{D*j ^ DK) e{D*j ^ D*K) 



(5) 



where the N{D*j) are the signal yields and the e{Dlj) 
are the corresponding efficiencies, and 1=1,4. We note 
that the only difference between numerator and denomi- 
nator final states is the presence of an extra 7r° from the 
D* decay. 

Assuming a constant total width, the yields are ob- 
tained by fitting the DK and D*K mass spectra using 
the same i:'*i(2710)+ and D*j{2860)+ parameters, and 
are summarized in Table IV. Efficiencies are evaluated 
using Monte Carlo simulations, and only the ratio of ef- 
ficiencies 



€r = 



e{Dlj ^ D*K) 
e{D*j ^ DK) 



(6) 



is involved in the measurement. This ratio is consistent 
with being uniform as a function of DK{D*K) mass and 
its values are reported in Table IV. 

Systematic uncertainties are summarized in Table V. 
The D*j parameters have been varied within their statis- 
tical and systematic errors, the p* cut has been changed 
to 3.1 and 3.5 GcV/c, the cos 9k cuts have been increased 
to —0.6. The systematic error arising from Monte Carlo 
statistics appears as the error on the ratio between the 
efficiencies. The error on the D* branching fractions is 
obtained from Ref. [14]. The shape of the background 
has been modified by adding an extra term in the ex- 
ponential; its contribution to the total error is found to 
be negligible. Finally the D* Am signal region has been 
reduced to ±2a. 

We obtain the following ratios of branching fractions: 



n = 



r-2 



where D*° 



r3 



B{Dti{2710)+ ^ D*°K+) _ 
B{D*^i2710)+ ^ DOK+) ~ (7) 

0.88±0.14stat±0.14syst 

i3(D*;(2860)+ ^ D*"K+) _ 
B{D;j{2860)+ ^ DOK+) ^ (8) 

1.04±0.17stat±0.20,y,t 
D°n°, and: 

g(D;i(2710)+ D*+K'^g) _ 
B{Dt,{2710)+^D+Kl) ~ (9) 

1.14±0.39stat±0.23syst 



BiDlj{2860)+ ^ D*+K^s) 
B{D*j{2860)+ D+K%) ~ 
1.38±0.35stat±0.49syst 



(10) 



where D*+ D+tt^. 

Averaging ri, and r2, we obtain: 



B{D*^{271Q)+ ^ D*K) 
B(£>*i(2710)+ ^ DK) 



0.91±0.13stat±0.12syst 
(11) 



B{Dlj{28m)+ ^ D*K) 

B(I?:,(2860)+ - DK) = '-'^ ^ °-'^«tat ± O.lQsysf 

(12) 

We also make a test of isospin conservation. We use the 
yields obtained from the fit to the appropriate mass spec- 
tra and correct for efficiency and branching fractions. We 
obtain, within the errors, similar rates for resonance de- 
cays to D°K+ and D+K° as well as for decays to D*^K+ 
and £)*+/r°, as expected from isospin conservation. 

We now compare these results with theoretical expec- 
tations. In the work of Ref. [15], for = 1~, two 
different quark model assignments are proposed for the 
D|i(2710)+: the I = 2 groimd state, l^Di, and the Z = 
ffist radial excitation, 2''5i. In the first case the ratio 
B(D*i(2710) D*K)/B{DI^{271{)) -> DK) is expected 
to be 0.043 ±0.002, in the second case 0.91 ±0.04. In this 
framework the D*]^(2710)+ can be identified as the first 
radial excitation of the D*(2112). The same assignment 
is derived from Ref. [6] where a mass of 2711 MeV/c^ is 
predicted for the state 2^51. However, in this case the 
expected ratio is 3.55, in significant disagreement with 
the measured value. 

In Ref. [16], in the framework of chiral doublers, 
= 1~ states are expected at masses of 2632 and 2720 
MeV/c2. 

The observation of £)*^(2860)+ D*K rules out the 
— 0+ assignment suggested by Refs. [5, 6]. In Ref. [7] 
the — 3^ assignment is proposed; however the pre- 
dicted B{Dlj{28m) D*K)/B{D*j{2860) DK) is 
0.39, which differs from our measurement at the level 
of three standard deviations. A better agreement is ob- 
tained if we compare with the calculations from Ref. [17] 
which expects a ratio of 0.6. 

As to the possible interpretation of the £>sj(3040)+ 
state, we note that Ref. [18] predicts two = 1+ radial 
excitations at 3082 and 3094 MeV/c^. 



IX. CONCLUSIONS 

In summary, in 470 fb"""^ of data collected by the BaBAR 
experiment, we observe the decays of the D*i{2710)^ and 
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TABLE IV: Information related to the evaluation of the ratio of branching fractions for the D*j resonances. 



Decay 




K 




B.F.{%) 




D:i(2710) + 
D:i(2710) + 
J (2860)+ 
i?: J (2860)+ 


^ D"K+ 
^ D-'"K+ 
^ D"K+ 
^ D*''K+ 


6469 ± 425 
1247 ± 173 
1826 ± 158 
415 ± 55 


0.353 ± 0.005 


61.9 ± 2.9 


0.88 ± 0.14 
1.04 ± 0.17 


D:i(2710) + 
13:1(2710) + 
£1:^(2860)+ 
D:j(2860)+ 


^D+K% 
^ D*+Kl 
^D+K% 
D*+K°s 


2442 ± 179 
258 ± 75 
781 ± 83 
100 ± 23 


0.301 ± 0.009 


30.7 ± 0.5 


1.14 ± 0.39 
1.38 ± 0.35 



TABLE V: Systematic uncertainties in the evaluation of the 
ratio of branching fractions. 



Ratio 


Dsj 


MC 


D* 


P* 


cos Ok 


Am 


Total 




parameters statistics 


B.F. 


cut 


cut 


cut 




ri 


0.030 


0.013 


0.042 


0.018 


0.120 


0.044 


0.14 


r2 


0.050 


0.015 


0.050 0.066 


0.175 


0.004 


0.20 


ra 


0.077 


0.033 


0.018 


0.042 


0.128 


0.160 


0.23 


r4 


0.106 


0.040 


0.022 


0.009 


0.328 


0.345 


0.49 



D*j(2860)+ to D*K and measure their branching frac- 
tions relative to DK. A new, broad D^j state is observed 
in the D* K mass spectrum at a mass near 3040 GeV/c^. 
Possible spin-parity assignments for these states are dis- 
cussed. 
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